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ABSTRACT 

We explore the effect of contamination of intermediate redshift Type la supernova samples by 
Type Ibc supernovae. Simulating observed samples of la and mixed Ibc/Ia populations at a range of 
rcdshifts for an underlying cosmological concordance model (fi m = 0.27, J7a = 0.73), we find that even 
small contamination levels, 2 — 5% may bias the derived J7a and f2 m towards larger values. We thus 
emphasize the need for clean samples of Type la SNe for accurate measurements of the cosmological 
parameters. We also simulate a SN sample similar to the fiducial SNAP detected distribution (Kim 
et al. 2004), but include Ibc contamination. For this distribution we illustrate the effect of Ibc 
contamination on the distance modulus vs. redshift diagram for low and high precision measurements. 
Subject headings: cosmological parameters — cosmology: observations — supernovae: general - 
surveys 



1. INTRODUCTION 

Type la supernovae (SNe) are considered astronom- 
ical standard candles and have been used to measure 
the geometry and dynamics of the Universe IjRiess et al.l 
119981 iPerlmutter et alJll99Hfl . This exceedingly difficult 
measurement relies on observations of local supernovae 
samples to calibrate relationships between absolute peak 
magnitude and more easily obtai ned observab l es suc h as 
light-curve width, color, etc. (sec Lcibundeut ( 2001) for 
a review). These calibrations are then applied to the 
observables in higher redshift samples to derive distance 
moduli. Due to the extreme faintness of intermediate to 
high redshift SN and the shifting of identifying spectral 
features outside the optical window, spectroscopic clas- 
sification becomes difficult, and photometry must often 
suffice. Fortunately, there is a clear observational differ- 
ence between hydrogen-rich (Type II) and hydrogen-poor 
(la, lb, and Ic) SNe, namely, the UV deficit (UV radia- 
tion is blocked in hydrogen-poor SN explosions by metal 
lines) . This property can be used to distinguish between 
Type I and II SNe explosions. To date there is no photo- 
metric technique to distinguish between Type la and Ibc 
SN (although see Gal- Yam et al., 2004), and one must 
perform a spectroscopic classification. The characteris- 
tic spectroscopic difference between Type la and Ibc SNe 
is the deep absorption trough at 6 140 A in la spectra, 
which is due to the blueshifted Si II AA6347, 637lA fea- 
ture. Unfortunately this is redshifted out of the optical 
passband for z > 0.5, making classification of Type I SNe 
difficult. Alternative spectral discriminants between la 
and Ibc SNe are Si II A41 30 (e.g. Coil et al . 2000), Fe II 
A4555 and Mg II A4481 l|Barris et al.ll200l . Spectro- 
scopic identifications are made using template-matching, 
but no statistical informa tion is given about misidentifi- 
cations with this method ijTonrv et al.ll2003t iRiess et al.l 
2004b] iBarris et al1l2004fl . 

Assuming the same relative la and Ibc rates in nearby 
SNe samples, Riess et. al (2004) and Barris et al. (2004) 
reasoned that the rate of contamination from Type Ibc 
SN in high redshift samples should be less than ~ 10%. 
However, the star formation rate increases rapidly at up 
to z ~ 1 — 2, and this is likely to have a significant effect 



on the relative rates of SN types. As Type Ibc SNe come 
from massive stars (M > 20 — 30 M Q ), their rate depends 
on the star formation rate in the previous 10 Myr. It is 
well established that the global star formation rate in- 
creases at higher redshifts. Hence, assuming a delay for 
la explosions, we expect that the rate of Type Ibc SNe 
increases over the short term, and that the contribution 
from Ibc SNe becomes more important at higher red- 
shift. Motivated by this, we investigate here what effect 
contamination by Type Ibc SNe has on cosmological pa- 
rameters derived from la SNe samples. The optimal way 
to determine the systematic effect in a particular survey 
is through simulations using the diagnostic machinery as- 
sociated with said survey. However, the general nature of 
the effect can be investigated robustly with simulations 
such as the one we present here. 

2. THE GLOBAL STAR FORMATION RATE 

The global star formation rate peaked some- 
where between z ~ 1 — 2 a nd has be e n de- 
clining to the present da y JLillv et al.l [T996t 
Madau. Pozzetti. fe DickinsoiJ 119981: iSteidel et all 
J^9^nTTomDS on^^fevmann. fc Storrie-Lombardin27)0lT 
Bouwens et al.ll2003D . Between the local Universe and 
z« 1.5, the star formation rate increased to 10 — 20 
times the rate in the current epoch. This significant 
increase in star formation rate should have a substantial 
effect on the relative rates of SN types. If the rates 
are about 1 Type Ibc SN for every 9 of Type la, when 
the star formation rate was 10 times its present value, 
the rates of Type Ibc and la SN should be about 
equal; for higher star formation rates, Type Ibc SN will 
outnumber those of Type la. Here we assume that the 
rate of Type Ibc SN explosions does not evolve with 
metallicity, therefore it depends only on the number 
of massive stars formed in the last 1 — 10 Myr, and 
thus scales linearly with the recent star formation rate. 
We also assume that the rate of Type la SN should 
be relatively constant after the "time delay" until the 
first explosion, which is somewhat speculative. It is 
possible that the rate of Type la SNe correlates with 
the star formation rate to some degree, but in any case 
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Fig. 1. — Magnitude distribution function for Type la and Ibc 
SNe assuming equal numbers of each and a single gaussian distri- 
bution for Type Ibc SNe. 
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Fig. 2. — Magnitude distribution function for Type la and Ibc SN 
assuming equal numbers of each and a double gaussian distribution 
with "bright" and "normal" Type Ibc SN. 



to a lesser degree than the core-collapse Type Ibc SNe. 
Whatever the actual degree, an increase in the global 
star formation rate will affect the relative ratio of Type 
la to Ibc SNe observed at that epoch, which leads us to 
the present investigation. Important constraints on the 
time delay for SNe la explosion are now b eing extracted 
from intermediate redshi ft SNe samples (Dahlc n~et al.l 
l2003 IStroker et al.ll200l . However, it is important to 
note that the general conclusions presented here do not 
depend on the Type la vs. Ibc rate with star formation 
rate, only that an increasing Ibc/Ia rate with redshift 
provided the motivation for this study. 

3. ESTIMATING THE EFFECT OF Ibc CONTAMINATION 

3.1. Absolute Magnitude Distributions 

We use the observed di stribution of uncorrected abso- 
lute magnitudes, from the lRichardson et alJ l)2002j) study 
of the Asiago Supernova Catalog. Pertinent here are the 
111 Type la SN, 5 luminous Type Ibc SN, and 13 "nor- 
mal" Type Ibc SN in the catalog. We will use their sin- 
gle gaussian fit to the uncorrected absolute magnitudes 



of the 111 Type la SN (M B = -19.16, a = 0.76) and 
their single gaussian fit to the Type Ibc SNe (Mb = 
-18.92, cr = 1.29); these are shown in Figure □ The 
number of Type Ibc SNe in the sample is small enough 
that a two-component model for the distribution of mag- 
nitudes is equally valid, with Mbi = —17.23, a\ = 0.62; 
M B 2 = -19.72,0-2 = 0.54, and a weight= 0.28. We run 
models with both the one and two gaussian distributions 
for Ibc magnitudes. 

3.2. Method 

We use a Monte Carlo simulation that creates a ran- 
dom redshift array from z = 0.5 — 3. Absolute mag- 
nitudes are assigned from the Type la and Ibc abso- 
lute magnitude distributions (described in the preced- 
ing section). We then include "observational scatter" by 
adding or subtracting flux from each object according to 
a gaussian distribution with a = 0.35 magnitudes. This 
is slightly larger than the m e an er ror of SNe between 
z = 0.5 - 2 in the lTonrv et~aTl l|2003ft sample (32 objects, 
mean error= 0.30). 

We "observe" objects by calculating the luminosity dis- 
tance for each object at its individual redshift, and "de- 
tecting" those with absolute magnitudes brighter than 
the detection limit, which we take as tub = 25. This 
should be roughly correct for redshifts up to z = 1. Be- 
yond this, rest-frame B corresponds most closely to J, 
and Mj = 24 is a more realistic limit. However, includ- 
ing this effect does not significantly affect the results in 
what follows. 

The Type la SNe that are "detected" are re-sampled 
to a Gaussian distribution with a center at what we 
take as the corrected Type la SNe absolute magnitude 
(-19.41; Leibundgut 2001), but with the dispersion equal 
to quadrature addition of the corrected dispersion of a 
typical low-redshift SNe sample (0.14; Barris et al. 2004) 
and an additional observational scatter. This simulates 
the corrections made to the Type la absolute magnitudes 
based on light curve width and color, but with a larger 
scatter to simulate observational error. Distance mod- 
uli are calculated for each object by assuming each has 
the absolute magnitude of a corrected Type la SN. Then 
the difference between these distance moduli and the dis- 
tance moduli for an empty universe is calculated. 

High redshift supernovae surveys utilize an analy- 
sis method based on the 'Multiwa velength' or 'Mul- 
ticolor Light-curve Shape Method' l|Riess et alJ l2004bt 
iBarris et alJ [2004) . This method uses a set of la 
lightcurves as templates for training, then concurrently 
fits the distance modulus, extinction, typically expressed 
in units of Ay, and A, a parameter which describes the 
lightcurve shape and corresponds to the difference in ab- 
solute magnitude of the SN and a fiducial SN la. 

The dominant correction comes from the correlation 
between light curve shape and peak magnitude, in the 
sense that fainter la SNe decline in luminosity faster than 
brighter ones. In order to discuss this quantitatively, 
we consider here the decline parameter, Amis, which is 
the difference between the peak magnitude and the mag- 
nitude at 15 days after maximum (e.g. Phillips 1993). 
The decline parameter varies with the observed restframe 
wavelength, and depends o n K-corrections, extinction, 
and r edshift (time dilation) ijNueent. Kim, fc Perlmutterl 
2002). The goal here is not to reproduce a particular 
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Fig. 3. — An example of detected SN redshift distributions (1000 input, 320 — 360 detected) SNe and varying the "observational scatter" 
from (a) 0.05, where one is dominated by la magnitude corrections, a = 0.14 and (b) 0.35, where the observational scatter dominates. 
Input (n M ,n A ) = (0.27,0.73). 



survey and its associated completeness and bias, but to 
investigate any possible systematic effect a contaminated 
SNe sample has on the derived cosmological parameters. 
Therefore, in this study we assume that the corrections 
(K, extinction, and time dilation) are accurately per- 
formed and the errors associated with these corrections 
are contained in the observational scatter parameter. 

To estimate typical m odifications to Ibc m agnitudes, 
we use the relation from lHamuv et ail (^996) for the B- 
band, M M ax = a + b[Am 15 (B) - 1.1], where b = 0.784, 
and we set a to the value we use for corrected peak la 
magnitudes, -19.41. We stress that the particular val- 
ues for MmaXi dispersions, etc. that we use are not 
significant, only that we consistently use the same values 
throughout the simulation. Then we need typical decline 
parameters for Ibc SNe. Most published lightcurves for 
lb and Ic SNe are not well-sampled in the period ±30 
days from maximum, but focus on the long-term behav- 
ior of the lightcurve. We have collected a sample of 5 
Ibc lightcurves with measurable Amis; these are sum- 
marized in Table 13.31 We also note that we h ave not 
listed the lb SN 1983N llClocchiatti et alJ 1199^1 and Ic 
SN 1983V ijClocchiatti et alJll997j) that have been shown 
to have similar lightcurves to SN 1993J, although the 
number of photometric points in the first ±30 days is 
small. 

3.3. Decline Parameters 

Decline parameters for these 5 Ibc SNe span a wide 
range of values. In light of this, we select reasonable 
ranges for the decline p arameter consi d ering the sample 
of la SNe presented in lHamuv et alJ 1^)96). We em- 
phasize that this is valid even if decline parameters for 
Ibc SNe are distributed differently than la decline pa- 
rameters, as Ibc SNe that are not weeded out will be 
fit to la templates. We explore four treatments assum- 
ing: (1) no modifications, (2) all Ibc SNe are matched 
to templates with small decline parameters (mean=0.9, 
dispersion=0.1), (3) all Ibc SNe are matched to tem- 
plates with large decline parameters (mean=1.5, disper- 
sion=0.3), and (4) Ibc SNe are matched with decline pa- 
rameters with a mean of 1.1 and a gaussian dispersion of 
0.3. This leads to mean assumed absolute magnitudes of 
(1) -19.41 (2) -19.57, (3) -19.10, and (4) -19.41 (same 



SN Type AA/i 5 M B , max B max - V max Rcf 



1990B Ic 0.964 (B) - (1), (2) 
1.14 (R) 

1992ar Ic -19.7 ±0.5 0.4 ±0.2 (3) 

19941 Ic 2.07 ± 0.03 (B) -17.68 ±0.73 0.15* (4) 

1993J lib 2.5(B) -17.3 ±0.6 0.18* (6) 
0.95 (V) (6) 
0.8 (R) (6) 

1998dt lb 0.2-0.3 (R) - - (2) 

1999di lb 0.4 (R) - - (2) 

1999dn lb 0.5 (R) - - (2) 

1999ex Ic -17.4 0.42 (7) 

1999bw Ic 1.0 (B) -18.88 OA (8) 



""Corrected for extinction, (1) Clocchiatti et al. (2001), (2) Math- 
eson et al. (2001), (3) Clocchiatti et al. (2000), (4) Clocchiatti et 
al. (1996) (5) Richmond et al. (1996), (6) Richmond et al. (1994), 
(7) Stritzinger et al. (2002) (8) Galama et al. (1998) 



as (1)). 

We compare these simulated data to a grid of mod- 
els with a range of fi m and fl\ values, and calculate 
the corresponding grid of x 2 values. Probability curves 
ar e calculated from t his x 2 grid following the method 
of iTonrv et alJ l)2003|) . That is, the minimum in the x 2 
array is found, and a grid of probabilities is calculated 
where prob=exp(— 0.5 x (x 2 — min(x 2 ))- The levels of 
68%, 95%, and 99.5% are found by integrating (sum- 
ming) the probability array and identifying where 68%, 
95%, and 99.5% of the total probability are located. This 
is only valid when the minimum reduced y 2 is close to 
1 (as it is in the SNe sample analyzed in ITonrv et alJ 

4. RESULTS 

4.1. Malmquist Bias with Ibc SNe 

In Figure 01 we show the change in the measured mean 
distance modulus (compared to an empty universe or 
eternally coasting model) from our simulation for input 
Ibc/Ia ratios of 0.1, 0.5, and 1.0 and an underlying con- 
cordance model with (Q m ,f2A.) = (0.27,0.73). An inter- 
esting effect can be seen, which is solely due to the obser- 
vational apparent magnitude limit. At redshifts less than 
~ 0.8, the mean distance modulus is fainter than the pure 
la case, and at redshifts greater than this it is brighter. 
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Fig. 4. — Change in the mean measured distance modulus for 
Ibc/Ia ratios of 0.1, 0.5, and 1.0. This is shown only for purposes 
of illustration. 



This is because of the luminosity distribution of Type Ibc 
SNe, with the majority being fainter than Type la SNe. 
At lower redshifts, Ibc contamination pushes the mean 
fainter, while at higher redshifts, these fainter SNe are 
below the detection limit, and only the brightest SNe 
arc detected, pushing the mean to apparently brighter 
magnitudes. The transition point is set by the apparent 
magnitude limit, here taken as tub = 25., and moves to 
higher redshift at fainter apparent magnitude limits, and 
vice versa. This is shown only for illustrative purposes, 
and no conclusions are derived from this diagram. 

4.2. Stochasticity 

The number of detected SNe affects the spacing of the 
probability curves, but more importantly, it affects the 
variation in the location of the \ 2 minimum from run to 
run. As expected, with smaller numbers of detected SNe, 
the curves are unstable, and "jump around" from run to 
run. As the number of input SNe is increased, and thus 
the number of detected SNe also increases, the position 
of the curves is more stable. 

As an example of this, for 100 trials using 500 Type la 
SNe uniformly and randomly distributed between z= 
0.5 — 2.0, with a mean of 190.04 detected SNe, the mean 
f2 m , Q\ values at the location of the y 2 minima were 
0.269 and 0.733, with standard deviations of 0.103 and 
0.181. The mean minimum reduced x 2 was 1.01, with a 
standard deviation of 0.11. Inputing 5000 Type la SNe in 
100 trials with a mean of 1884.07 detected SNe decreased 
the standard deviations to 0.034 and 0.064, with mean 
n mi Q\ values of 0.27 and 0.73. The mean minimum 
X 2 /N was 1.007, with a standard deviation of 0.031. 

We investigated what effect Type Ibc SNe would have 
on the derived probability curves by assuming values of 
0.27 and 0.73 for Sl m and f^A (the 'cosmological concor- 
dance' model) and distributing 2000 Type la SNe uni- 
formly and randomly between z = 0.5 — 2. In one run 
of the simulation, 724 of the 2000 were "detected", and 
the probability curves corresponding to the y 2 analysis 
are shown in Figure Efa). Using these same 2000/724 
Type la SNe, we first added 100 Type Ibc SNe with 
magnitudes distributed randomly according to the uni- 
modal distribution presented in § 13.11 11 of these were 
"detected" , and the result of this simulation is shown in 



Figure Efb) . In Figure OJc) we show the result of us- 
ing 200 Type Ibc SNe (input ratio of 10%), of which 26 
were detected. The reduced x 2 increased from 0.97 for 
the pure la sample, to 1.13 for the sample used in Fig- 
ure E{b), to 1.36 at a contamination level of only 5%, 
shown in Figure p^c). 

In these simulations, relatively low contamination rates 
of 26/742=4% significantly affect the location of the 
probability curves in the £Im^a plane, moving both 
f2M,^A to higher values. However, as the fraction of 
detected Type Ibc SNe increases, the minimum reduced 
X 2 value also increases. At some point, this minimum 
reduced x 2 value should be large enough that one would 
correctly conclude that the model is not a good fit, and 
the conversion from the x 2 grid to probability curves 
would not be, and should not be, attempted without 
first renormalizing the errors. Then the probability con- 
tours become larger, and should encompass the correct 
values. However, in this renormalization it is im- 
plicit that the model is a good fit, and the random errors 
are underestimated. In this example, this is incorrect. 

The spacing of the probability curves changes with the 
value of the observational scatter and the total number 
of SNe detected. But the stability of the location of the 
center of the curves depends on the number of detected 
SNe, varying stochastically. We find that for large sam- 
ples with relatively accurate measurements, even a small 
contamination (2 — 3%) of Type Ibc SNe can significantly 
bias the result. In 100 trials with an input of 5000 SNe 
uniformly and randomly distributed between z= 0.5—2.0, 
an observational scatter of 0.35, an input Ibc/Ia ratio of 
0.1, and a bimodal magnitude distribution (resulting in 
a mean of 1891.05 Type la and 54.61 Type Ibc SNe de- 
tected), the mean f2 m ,fiA values for the location of the 
X 2 minimum were 0.334 and 0.826, with standard devia- 
tions of 0.038 and 0.069. The minimum reduced \ 2 was 
1.156, with a standard deviation of 0.048. In this case, 
the underlying values of f2 m = 0.23, = 0.73 are out- 
side the la curves, but inside the 2a curves, 68% of the 
time. 

It is worth noting that the minimum y 2 is at 
(n m ,fi A ) ~ (1.3,0.6) for the lTonrv et al. ( 2003§ sample, 
before applying external constraints. In fact, without 
external constraints, the "concordan ce model" values are 
outs ide the la co n tours. However, in lRiess "elaTl (I2004bf) 
and iBarris et alJ (|2004) , the error contours are signifi- 
cantly larger, and encompass the concordance model val- 
ues for the density parameters. 

4.3. Unimodality vs. Bimodality 

We also compared the effects of single and double gaus- 
sian Ibc magnitude distributions. These results are sum- 
marized in Tabled With a single gaussian distribution, 
stochastic effects are more important, as illustrated by 
the larger standard deviations for the mean recovered 
fiju, ^A, and minimum x 2 /N. In terms of systematic 
error in the recovered cosmological parameters, if the 
double gaussian distribution for Ibc magnitudes is more 
accurate in describing the underlying distribution, con- 
tamination by Ibc SNe will bias the results less severely 
than the single gaussian case. 

4.4. Effect of Decline Parameters 
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Fig. 5. — Input of 2000 la uniformly and randomly distributed from z= 0.5 — 2 and a single gaussian distribution of (b) 100 Ibc and 
(c) 200 Ibc. (a) 724 detected Type la SNe. The minimum reduced \ 2 is 1-03. (b) 724 detected Type la SNe, 11 Ibc SNe. The minimum 
reduced x 2 is 1.13. (c) 724 detected Type la SNe, 26 detected Ibc SNe. The minimum reduced x 2 is 1.36. 
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Fig. 6. — See text for details. From left to right: (a) 231 la SNe, (b) 231 la, 15 Ibc with no modification to the Ibc magnitudes, (c) 
231 la, 15 Ibc with modification with small decline parameters, (d) 231 la, 15 Ibc with modification with large decline parameters, (e) 
231 la, 15 Ibc with modification with medium decline parameters. There is no significant difference between scenarios with this range of 
Ibc magnitude modifications. 



TABLE 1 

A COMPARISON OF THE MAGNITUDE DISTRIBUTIONS FOR IBC 
SUPERNOVAE. RESULTS FROM 100 TRIALS. 



single gaussian double gaussian 



negligible compared to the magnitude difference between 
the detected la and Ibc SNe. 

4.5. Colors 

From the sample of SNe compiled in Table EH we find 
that the B max — V max colors of Ibc SNe have a tendency 
to have redder colors than la SNe, which typically have 
B m ax — V max ~ 1- However, the slope of the relation 
of B — V color with time f rom maximum is simila r, as 
illustrated in Figure 15 of Stritzin ger et al.l l)2002|) for 
SN1999ee (la) and SN1999ex (Ib/c). If the maximum 
was not observed, a shift of ~ 15 days would place the 
two B — V curves on top of one another. 

If the B max — V m ax color is correctly measured through 
an accurate determination of the date at maximum, then 
an interloping Ibc SN will appear redder than the average 
la SN. If this Ibc SN has a large decline parameter, it will 
be "expected" to be fainter, and fainter la SNe are redder 
at maximum (see Figure 1 of Riess, Press, & Kirshncr 
(1996) for examples). However, if this Ibc SN has a small 
decline parameter, it will be "expected" to be brighter, 
and brighter SNe are bluer. Therefore, the red color will 
be interpreted as due to extinction. 

As an example, let us consider SN 1998bw 
ijGalama et al.l 11998^1 . the first confirmed super- 
nova/ GRB. This supernova had a decline parameter 
Am 15 = 1.0, Bmax - Vmax of 0.4 magnitudes, and 
Ms maT = —18.88. If the intrinsic color was assumed to 



mean n a/ 0.409 ±0.126 0.344 ± 0.094 

mean Q A 0.974 ± 0.243 0.839 ±0.176 

mean minimum x 2 /N 1.393 ±0.242 1.225 ±0.110 

mean number of SNe la detected 228.77 226.78 

mean number of SNe Ibc detected 12.46 10.26 



We explore four treatments assuming: (1) no modifi- 
cations to Ibc magnitudes, (2) all Ibc SNe are matched 
to templates with small decline parameters (mean=0.9, 
dispersion=0.1), (3) all Ibc SNe are matched to tem- 
plates with large decline parameters (mean=1.5, disper- 
sion=0.3), and (4) Ibc SNe are matched with decline pa- 
rameters with a mean of 1.1 and a gaussian dispersion of 
0.3. This leads to mean assumed absolute magnitudes of 
(1) -19.41 (2) -19.57, (3) -19.10, and (4) -19.41 (same 
as (1)). 

In Figure we compare the results for a single case for 
these four treatments. The number of detected SNe was 
231 la and 15 Ibc. Assuming Ibc SNe are treated as la 
SNe or the magnitudes are not modified does not have a 
significant effect. The corrections to the magnitudes are 
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Fig. 7. — Including a changing star formation rate with redshift 
with a single gaussian Ibc magnitude distribution. In (a), the rate 
of Type Ibc to la SNe follows the dependence 0.05(1 + z), while 
in (b) it varies as 0.1(1 + z), the same redshift dependence but 
with a different normalization. The number of detected la and Ibc 
SNe was (a) Ia= 318, Ibc= 20, and (b) Ia= 287, Ibc= 41. The 
minimum x 2 /N is (a) 1.57 and (b) 2.27. 



Fig. 8. — Including a changing star formation rate with redshift 
with a two gaussian Ibc magnitude distribution. The la distribu- 
tion is identical to that in Figure [7] In (a), the rate of Type Ibc 
to la SNe follows the dependence 0.05(1 + z), while in (b) it varies 
as 0.1(1 + z), the same redshift dependence but with a different 
normalization. The number of detected la and Ibc SNe was (a) 
Ia= 318, Ibc= 15, and (b) Ia= 287, Ibc= 32. The minimum x 2 /N 
is (a) 1.21 and (b) 1.52. 



be B max - V m ax = 0, then E(B - V) = 0.4, A v = 1.24, 
and Ab = 1.64. The observed apparent rest-frame 
B magnitude would be corrected by approximately 1.6 
magnitudes, and the derived distance modulus would be 
too large by more than 1 magnitude. Put another way, 
this object would appear 1 magnitude brighter than an 
accurately corrected la SN. 

This effect should act in a random way, unless the 
peak magnitudes, colors, and light curve shapes of Ibc 
SNe are correlated. With the current available sample in 
the literature, we are unable to determine the relation- 
ship between these parameters. However, we note that 
some type la SNe have B max — V max colors similar to 
those of Ibc SNe, for example, 1990Y, 1992K, and 1993H 
ijHamuv et alj fl996). If these or similar lightcurves are 
used as templates in the fitting procedure, then the ar- 
gument in the previous paragraphs is negated. For our 
example Ic SN 1998bw, the "corrected" mag would be ap- 
proximately 0.6 magnitudes too faint. In this case, the 
effect of intrinsic B max — V max colors, which do appear 
to be different for la and Ibc SNe, can only be correctly 
assessed by testing the template fitting codes of individ- 
ual surveys. We urge investigators to consider testing the 
robustness of their photometric and spectroscopic fitting 
to Ibc interlopers. 

4.6. SFR(z), Ibc/Ia(z) 

We also investigated a Ibc/Ia ratio that changes with 
redshift. In Figure we present examples of including 
a dependence of the Ibc/Ia ratio with redshift, mimick- 
ing a dependence on the star formation rate. Due to 
the severe effect of Ibc contamination on the probability 
curves even for small levels of contamination, we were 
forced to keep the relationship shallow (1 + z) and with 
a small Ibc/Ia normalization ratio (< 0.1). In Figure [7| 
the minimum reduced x 2 was 1-5- However, for the case 
shown in Figure [H] the minimum reduced \ 2 was 2.6, 
and the simple calculation of probability curves we have 
shown is not valid. A more sophisticated treatment us- 
ing the estimated la and non-la SNe rates with redshift 
from Dahlen et al. (2004) could be performed, but the 
measurements are so sensitive to Ibc contamination that 
a more complex treatment is unwarranted. 



4.7. Comparison to the SNAP Redshift Distribution 

iKim et all ((2004) reviewed systematic uncertainties in 
SNe studies to determine the cosmological parameters. 
In their study, a fiducial number-redshift distribution for 
SNe was used, an observational scatter of 0.1 magnitudes, 
and a dispersion about the corrected la peak magnitudes 
of 0.1 magnitudes. We would like to use the same SNe- 
rcdshift distribution to see what effect Ibc contamination 
would have. However, including the redshift distribution 
of detected SNe is problematic, as instead of specifying 
detected SNe, we assign an input distribution, then "ob- 
serve" those above our apparent magnitude limit. To ac- 
curately determine the effects of Ibc contamination, an 
apparent magnitude limit must be included, thus some 
"observation" must be simulated. To match the SNAP 
fiducial distribution, we tailor the input distribution such 
that the "observed" distribution is sufficiently similar to 
the fiducial SNAP distribution. 

We attempt to reproduce the SNAP fiducial distribu- 
tion's general increase in numbers from z = 0.1 — 0.8, the 
peak near z=0.8, and the slope of the decline in numbers 
from z = 0.8-1.7, dlogN{z)/dz = -0.4. The solution of 
input to detected numbers is not general, but depends on 
the detection limit. For mg = 25., we find a reasonable 
agreement with an input of two power law distributions, 
700 SNe between z = 0.1 - 1.0 as N oc z 2 , and 10000 
between z = 0.8 — 1.7 as N cx z 5 . We then input 1000 
Ibc SNe between z = 0.8 — 1.7 as N tx z 5 , the same re- 
lationship as for the high redshift la population. Thus, 
we include an underlying ratio of 0.1 for Ibc/Ia SNe. An 
example of a final detected redshift distribution is shown 
in Figure EJa), an d the A(m — M) vs. z distribution 
in Figure OJb) . We emphasize that this is one run of a 
Monte Carlo code, and as in every random process, each 
individual answer is different. 

An important issue is clearly illustrated in Figure EH if 
the observational scatter and the dispersion about the 
corrected peak magnitudes is low, 0.1 magnitudes for 
each in this simulation, then Ibc SNe are easily distin- 
guished. Even when the number of SNe is low, the bright 
Ibc SNe that bias measurements of the cosmological pa- 
rameters are easily isolated. On the other hand, when 
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the observational scatter is large, large numbers of SNe 
are needed to identify outliers, in this case Ibc SNe. This 
is illustrated in Figure ITU1 

//such a low scatter and corrected dispersion is feasi- 
ble, then the need for spectra of SNe candidates to deter- 
mine Ia-ness is eliminated. Depending on the distribu- 
tion of SNe, such a survey could possibly rely instead on 
rejecting 3<r outliers from the mean or median distance 
modulus in a given redshift bin. Even if the observa- 
tional scatter is ~ 0.3, if the numbers are large, then 
er-clipping is still a valid option to remove outliers, how- 
ever, the constraints on the cosmological parameters will 
not improve, as the probability curves will not shrink. 
The important question becomes, is such precision fea- 
sible? The answer to this question is best left to survey 
designers. 

5. ON THE LIKELIHOOD OF Ibc CONTAMINATION 

As the progenitors of Type Ibc SNe are massive stars, 
we can expect that the majority have higher extinctions 
than la SNe and occur in gas-rich galaxies. Thus, ex- 
cluding SNe with high derived Ay values and those oc- 
curing in late-type host galaxies should reject more Ibc 
SNe than la and reduce contamination. However, the SN 
rate in different galaxy types as a function of redshift is 
not known. Although the local Ibc/Ia rate is ~ 0.1, we 
can expect this to increase to 0.5 — 1.0 around 2 = 1. 
Assuming no evolution in the luminosity functions, this 
means 1—3 bright Ibc SNe for every 10 la. For Ibc 
contamination to be a significant concern, 20 — 50% of 
these bright Ibc SNe must fail to be rejected through 
the possible combinations of light curve fitting, spectral- 
typing, host galaxy morphology, and extinction. The 
likelihood of this could be determined through detailed 
simulations using Ibc templates (e.g. 1999bw, 199 2ar) 
and rou tines such as ML CS2k2 (Ri ess et alJl2004q) . or 
BATM (|Tonrv et al.ll200l . 

6. CONCLUSIONS 

We find that even small contamination fractions (5%) 
may bias the measurement of VLm and f^A- However, we 
have also illustrated the ease with which Ibc SNe can be 
identified when the observational error and the dispersion 
about the corrected la peak magnitudes is low. Based on 
the calculations presented here, the straightforward effect 
of contamination from Type Ibc SN in intermediate-to- 
high redshift SN samples is to reduce the x 2 goodness-of- 
fit to the true cosmological parameters. This manifests 
itself as an increase in the minimum value of the reduced 
X 2 grid. For a flat universe with Qm — 0.27, Vt\ = 0.73 
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even small numbers (5%) of Type Ibc SNe may push both 
values upwards. Erroneously small "probability" curves 
can be generated by assuming the model is a good fit and 
that the minimum reduced \ 2 is close enough to 1 for a 
meaningful conversion to a probability grid. However, if 
the observational errors are large, it is difficult to distin- 
guish between a minimum reduced \ 2 that is larger due 
to a mismatch of the model and the observations due to 
contamination or from simple stochastic effects. 

If the observational errors are large enough (> 0.3), 
then a small contamination fraction (5%) by Ibc SNe can 
shift the location of the x 2 minimum without making x 2 
itself large. In such a case the error contours will be large, 
so while the contours are not strictly statistically correct, 
they are still meaningful in the sense that the underly- 
ing cosmological parameters have close to a 99% chance 
of being within the 2cr contours. If the observational 
error approaches the scatter in the corrected 'standard 
candle' la magnitudes, the errors bars will shrink. How- 
ever, with the same amount of Ibc contamination, as the 
observational error steadily decreases, the reduced x 2 be- 
comes steadily greater than 1. Then the conversion from 
a x grid to probabilities is not valid. One may assume 
that the random errors are underestimated and renor- 
malize; the error contours may then become large enough 
to encompass the underlying Hm^a values. But if this 
precision is achieved, outliers should already be readily 
discriminated in the A(m — M) vs. z plane, and a sim- 
ple 3 — 5er clipping routine may suffice. Simply having 
large numbers of SNe with the precision of current sur- 
veys would also allow one to use er-clipping to remove 
Ibc outliers, however, this will not achieve significantly 
greater precision in the measurement of cosmological pa- 
rameters. 
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Fig. 9. — Observational scatter of 0.1 magnitudes. Top panels: 188 detected la, 9 Ibc, bottom panels: 1922 detected la, 64 Ibc. Stars in 
the right panels indicate Ibc SNe. 
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Fig. 10. — Observational scatter of 0.35 magnitudes, top panels: 207 detected la, 10 Ibc, bottom panels: 2029 detected la, 80 Ibc. Stars 
in the right panels indicate Ibc SNe. 



